Nonequilibrium steady state fluctuations in actively cooled resonators 
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We analyze heat and work fluctuations in the gravitational wave detector AURIGA, modeled as 
a macroscopic electromechanical oscillator in contact with a thermostat and cooled by an active 
feedback system. The oscillator is driven to a steady state by the feedback cooling, equivalent to 
a viscous force. The experimentally measured fluctuations are in agreement with our theoretical 
analysis based on a stochastically driven Langevin system. The asymmetry of the fluctuations of 
the absorbed heat characterizes the oscillator's nonequilibrium steady state and reveals the extent 
to which a feedback cooled system departs from equilibrium in a statistical mechanics perspective. 



Cold damping feedback efficiently reduces the ther- 
mal noise motion of an oscillator by applying a viscous 
force. Since its first application in electrometers [J], it 
succeeded in a vifide variety of devices 0, 01 , from nano 
to macroscopic resonators, and in a variety of implemen- 
tations, including both optical and electrical forces. In 
basic research, the cold damping is considered in order to 
reduce below the level of intrinsic quantum fluctuations 
the uncertainty due to thermal noise of the position of 
macroscopic bodies [3], and to improve the behavior of 
gravitational wave detectors In this work we exper- 
imentally investigate the fluctuations of thermodynamic 
quantities of a cold-damped electromechanical oscillator: 
the resonant-bar gravitational wave detector AURIGA 
In particular we verify that they are consistent with 
recent theories of nonequilibrium phenomena. 

After the seminal works of Ref. 0] , which introduced 
the Fluctuation Relation (FR) concerning the Probabil- 
ity Density Function (PDF) of the entropy production 
rate in nonequilibrium systems, a large number of pa- 
pers has been devoted to similar problems (see for in- 
stance Ref. for a review). One finds that the FR 
for some properly identified observable (called dissipation 
function) is quite generally valid in systems of physical 
interest [9] . After the experimental evidence obtained for 
dragged colloidal particles [lo|, electrical circuits 11] and 
mechanical oscillators [l^, the FR has become a stan- 
dard tool to characterize nonequilibrium systems. Here, 
following Ref. [13] we focus on the FR for the PDF of 
the power necessary to maintain a dissipative system in 
a nonequilibrium steady state (NESS). A specific FR, 
obtained for the fluctuations of the injected power in a 
stochastically driven Langevin system [ij], was recently 
confirmed in a simple electrical realization of that model 
[lit . Actually, this FR accurately fits the fluctuations of 
the injected power in wave turbulence as well [l6| . We 



show that also the AURIGA detector, which is main- 
tained in a NESS by an external driving in a feedback 
cooling scheme, can be described as a mechanical oscil- 
lator forced by a stochastic driving. We then analyze 
its behavior and demonstrate that: a) the statistics of 
its thermodynamic variables show a characteristic asym- 
metry between positive and negative fluctuations; b) the 
statistics of the injected power are in agreement with the 
FR of Ref. [14]. These results reveal the extent to which 
a feedback cooled system departs from equilibrium in a 
statistical mechanics perspective and prove the limits of 
usual assumption that cold-damped oscillators at tem- 
perature To are equivalent to higher-loss ones, in ther- 
modynamic equilibrium at a temperature Toff < Tq. 

AURIGA is based on a 2.2 x 10^ kg, 3 m long bar made 
of a low- loss aluminum alloy (A15056), cooled to liquid 
helium temperature Tq = (4.6 ±0.2) K. The fundamental 
longitudinal mode of the bar, sensitive to gravitational 
waves, has effective mass M=l.l x 10^ kg and resonance 
frequency ujo/2tt ~ 900 Hz. The bar resonator motion is 
detected by a capacitive transducer followed by a double 
stage dc-SQUID amplifier; the displacement sensitivity 
is about 5 x 10^^^ m/\/Hz over a ~ 100 Hz bandwidth 
around ujq, largely limited by thermal noise. The detec- 
tor can be modeled by three coupled low- loss resonators: 
two mechanical ones (the bar and a plate of the capaci- 
tive transducer) and an LC electrical one [13, [13 • Their 
dynamics is described by three normal modes at separate 
frequencies, 865.7, 914, 953 Hz, with quality factors re- 
spectively of 1.2 X 10^, 0.88 X 10^, 0.77 x 10^, determined 
by mechanical losses in the bar and the transducer and by 
dielectric losses in the electrical components. Each mode 
is modeled as a RLC series electrical oscillator with an 
effective inductance L, capacitance C and resistance R, 
which assume different values for the 3 oscillators (Fig.[T]). 
These modes are electromechanical rather than purely 
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FIG. 1: (Color online). The normal mode is approximated, 
around its resonance frequency, by a series- RLC circuit. The 
dc SQUID is represented as current amplifier. The observ- 
able is the current 7s, and the electronic feedback cooling 
is obtained by sending back a current Id which is a delayed 
copy of Is reduced by G <C 1. The SQUID output voltage is 
Vout = AIs with A = 2.6 lO'' Q.. 



mechanical, but each one collects a significant fraction of 
the energy of the two mechanical resonators. From the 
detector calibration we estimate that the energy injected 
by an impulsive excitation of the longitudinal mode of 
the bar is shared by the 3 modes in the ratio 48:36:16. 

To the sole purpose of improving the electronics sta- 
bility and easing the data analysis, AURIGA employs 
an electronic feedback cooling scheme on the detector 
(Fig. [T]), which is equivalent to a viscous force on the 
oscillators [19]. The dynamics of each electromechanical 
oscillator is described by the equations: 
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where the current Ig is the observable, q is the charge on 
the capacitor, / = ^^^^ the current through the induc- 
tance L, Vd is the voltage at the node where the feed- 
back is applied and Li„ is the input inductance of the 
SQUID amplifier. In thermodynamic equilibrium, each 
oscillator is driven by the stochastic voltage: Vxit) = 
y/2kBToR r{t) , where F is a Gaussian white process. 
This should hold even in our nonequilibrium case, since 
the feedback cooling concerns only 3 modes, out of the 
very many degrees of freedom of the thermal bath, and 
is not expected to significantly affect the thermal noise 
due to the bath. The noise due to the SQUID can be 
neglected at moderate feedback gains, as those used in 
this experiment. The very high quality factor of the os- 
cillators imply that the currents Is, Id and / oscillate 
at ujr — 1/VlC with amplitude and phase changing ap- 
preciably only on timescales of several cycles. Thus, the 
quasi-harmonic approximation I(t) = I(t) sm[u!rt + cl){t)], 
like the analogous ones for Is and Id, seems appropriate. 
Operatively this is implemented by considering the signal 
only in a narrow frequency band around ujr. Near reso- 



nance, a feedback force equivalent to a viscous damping 
can be obtained with the feedback current: 



Idit) GIs{t - td) 



(2) 



where td — ^j- and G <C 1. Equation |Tb| includes 
memory effects due to contributions from times t — td, 
because of the constraints (fTc)) and ([2]); in the quasi- 
harmonic approximation we have Is{t — td) — ujrQs{t)- 
Hence each oscillator obeys: 



L^^+Is{t) [R + Rd] 



dt 



c 



= V2kBToRr{t) (3) 



with Is{t) — ■ Here Rd = Guj^Lm expresses the 

viscous damping on the oscillator due to the feedback 
loop; the feedback efficiency is defined as g = Rd/R- 
The quasi-harmonic approximation is valid as long as the 
feedback damped oscillator is still a low loss one. In Eq. 
([3|) the driving is the same white process of Eq. H]): 
this is confirmed experimentally by the Lorentz-shaped 
jower spectrum of the current Is around the resonance 
19| . Equation ^ is not invariant under time reversal 
Ws = 1st I's — ~Is, t' — —t) and does not satisfy the 
Einstein relation. Nevertheless, it is formally identical 
to that describing an oscillator with damping R + Rd, 
in equilibrium at the fictitious "effective temperature" 
Toff = To/(l -I- 5). The discrepancy between Tgg and 
the thermal bath temperature Tq reveals the nonequilib- 
rium nature of the phenomenon. Hence, the feedback 
cooled oscillator is usually treated as an equilibrium sys- 
tem, with Teff derived from the experimental value of 
Is{t)'^ = 2^^jf^, even if no bath at T^s is present. 
Multiplying Eq. ^ by Is{t) and integrating be- 
tween t and t + T, in the quasi-harmonic approxima- 
tion we get an expression for the average power Pt — 
T It^^ -^s{t')VT{t')dt' injected by the stochastic thermal 
force during a time t: 

= AUr + {t')dt' (4) 



where AC/t- 



U{t + T)-U{t) 



U (t) being the stored energy: 



Uit)^hlUt) + \'^ = lLi^{t) (5) 

The term proportional to R represents the heat dissi- 
pated by the oscillator toward the bath while that pro- 
portional to Rd is the work done by the oscillator on the 
feedback: 



^ ^'^\d{t')Vd{t')dt' ii{t')dt' (6) 



Notice that the last identity is strictly valid only within 
the quasi-harmonic approximation, which relates both 
Is(t—td) and dls(t)/dt to the instantaneous current Is{t). 
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FIG. 2: (Color online). PDF (units of s/^ksT^a) of: (a) the 
time averaged energy difference A(7t, (b) the work Wt and 
(c) the heat Qr averaged at increasing values of the ratio 
r/Tefi- Data were collected by AURIGA in a 10 days time 
span. Dashed vertical line at 0.84 fcsTed/s corresponds to 
the mean value of Wt- The gray dotted lines are obtained 
by numerical simulation of Eq. Q for a 50 days time span. 
The discrepancies observed at short r between experimental 
and numerical data are within the uncertainty due to the 
experimental error in TeS- 

Further, if r = N'jj-, N integer, in the same approxima- 
tion we can also write: 

By energy conservation we obtain the heat Qr , absorbed 
by the oscillator from the bath and averaged in a time r: 

^ AUr + Wr (8) 

To study nonequihbrium properties, we focused on the 
lowest frequency mode out of the 3, which is well sep- 
arated in frequency from the other two and is thus our 
best approximation of a single oscillator. The sampled 
current Is (t) was processed via the AURIGA data anal- 
ysis and integrated over the resonance in a 10 Hz band- 
width to obtain the current amplitude Is{t) in the har- 
monic approximation. From dedicated calibration of AU- 
RIGA we measure L = (1.67 ± 0.01) 10""* H and Lin = 
(1.48 ± 0.01) X 10"'^ H. A first set of data covers a con- 
tinuous 10 days time span acquired in March 2008. They 



yield iVr/^ir = 865.7 Hz, T^s = (21.1±0.2) mK and decay 
time Tcff = (2.36±0.04) s; hence we estimate g = 207±10, 
R = (6.8 ± 0.5) X 10-^ r2 and G = (1.74 ± 0.06) x lO^^. 
The quality factor uJrTeff/2 ~ 6.5 x 10'^ is high enough to 
justify the quasi-harmonic approximation leading to Eq. 
^ . Figure [5^ and [^b show the PDF of the energy dif- 
ference AUr and of the work done by the oscillator WV 
averaged over growing times r: they are calculated from 
Eqs. ^ and ([7]) after dividing the experimental data in 
contiguous time intervals of duration r. Figure [2t shows 
the corresponding heat Qr exchanged by the oscillator 
with the bath averaged over the time r, computed via 
the energy conservation Eq. ([8]). The fluctuations of Qr 
are asymmetric, as expected for a NESS. Figure [2] shows 
also excellent agreement with numerical simulations of 
Eq. dSl). 

The PDF of AUr is symmetric with respect to zero 
as for an equilibrium oscillator. It has exponential tails 
which decay faster for longer r. The PDF of Wr is highly 
asymmetric. From Eq. ^ and (/^(t)) = 2kBT^a/L 
we infer that Wr is positive and has mean value ~ 
0.84 ksTcs/s. Hence, Qr takes negative values only 
for short integration times, with the characteristic time 
scale given by the cold damped oscillator decay time 
Toff — 2L/{R + Rd). For r ^ Tcff the contribution of 
the time averaged energy is negligible. So in the pres- 
ence of feedback [Rd > 0) there is a net heat transfer 
from bath to oscillator: this is the energy flux that feeds 
the NESS and makes the PDF of the heat asymmetric. 
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FIG. 3: (Color online). Comparison between theory (dashed 
lines) and experimental data at increasing values of the ratio 
t/tcS = 1.2, 2.0, 2.9, 3.7: the dataset are respectively 828, 
768, 729, 457 days long, a) Plot of second derivative D2(er); 
the experimental data of D2 for t/tcH = 2.0 are not shown for 
clarity, b) Plot of p(er). Vertical error bars on experimental 
points of both plots come from statistical uncertainty. The 
shaded areas on the theoretical curves represent the uncer- 
tainty due to the experimental errors on the parameters. Toff, 
Teff and To- 
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On the contrary, if the feedback were switched off, we 
would have Rd = and Wt = 0, hence Qr = by 
Eq. ([5]). In this case the PDF of Qr would be symmetric 
with respect to its (zero) mean value, as in Fig. [5^, but 
with Tefi = Tq. 

The PDF of the injected power is essentially iden- 
tical to that of Qt, shown in Fig. [2j;, since Pr « Qt 
when g ^ 1. Notice that of the two terms in Eq. ([4]), 
only MJt is responsible for the negative values of Pr- 
Thus, large positive values of Pt are dominated by the 
contribution of the dissipated power [the integral in ^] 
more than they are for small values of Pt . The transition 
between these two regimes affects the shape of the PDF, 
which has been calculated in Ref. 
integration times it obeys: 



lj| . In a limit of large 



f{er)= lini 



lnPDF(?^) 



-7(1 - 2~er) 



if e-r 

if Ir 



(9) 

where e-^ — Pt-L/ [kBT^R) is the reduced injected power 
and 7 = (i? + Rd)/L = 2/tcs- A remarkable singularity, 
located at = 1/3, is present in the second derivative of 
/(cr). In Fig. [3^ we plot the quantity £'2(er) = ^ g^i^^ 
evaluated from the output of AURIGA in the timespan 
May 2005/May 2008; here TeS = (22 ± 1) mK and Tcs = 
(2.4±0.2) s. A valley is clearly visible in the experimental 
data, which we interpret as precursor of the asymptotic 
singularity. The agreement with the asymptotic theory 
consistently improves as t/tcs grows. This indicates that 
the asymptotic relation Eq. ( |9]) h olds even in presence of 
a harmonic pinning potential [14| . 

Eq. ^ easily leads to the FR for the injected power, 
i.e. to the ratio between the probability of positive and 
negative fluctuations of e,-. If we define p(er) —r™aD 
T In , we have: 
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As shown in Fig. [31d, positive values of e-r are exponen- 
tially more probable than negative ones. Two conflicting 
features determine the details of the experimental curves: 
the agreement with the asymptotic theory improves with 
t/tcs, but the statistics blur at large values of ir since 
negative events are rarer. For this reason a slope change 
is clearly seen for small values of r/ TeS where the precur- 
sor of the singularity occurs as shown in Fig. [3^, while 
it becomes barely visible at longer integration times. 

In conclusion, we demonstrate that the actively cooled 
AURIGA detector is well described by the Langevin 
model of Eq. ([3]), which led us to evaluate the power 



Pr injected by the stochastic thermal force, the work WV 
done on the feedback and the heat Qr exchanged with 
the thermal bath. The statistics of Pr are consistent with 
Eq. ([9|) , and with the consequent nonlinearity of the FR. 
The fluctuations of Qr are asymmetric as expected for a 
NESS. 
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